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Abstract 

The Schrodinger equation with the PT— symmetric Hulthen potential is 
solved exactly by taking into account effect of the centrifugal barrier for any 
/-state. Eigenfunctions are obtained in terms of the Jacobi polynomials. The 
Nikiforov-Uvarov method is used in the computations. Our numerical results 
are in good agreement with the ones obtained before. 
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I. INTRODUCTION 

The Hulthen potential [1] is a short-range potential in physics. The potential form is 

— St 

V C\r) = -ZtS^^s, (1) 

where Z is a constant and 5 is the screening parameter. If the potential is used for atoms, 
the Z is identified with the atomic number. This potential is a special case of the Eckart 
potential [2], which has been widely used in several branches of physics and its bound states 
and scattering properties have been investigated by a variety of techniques [3]. 

The radial Schrodinger equation for the Hulthen potential can be solved analytically 
only for s— states (/ = 0) [1,4,5]. For I ^ 0, a number of methods have been used to find 
bound-state energies numerically [3,6-16] and analytically [17,18]. In this context, we present 
in this letter a Nikiforov-Uvarov method [19] within the frame of PT— symmetric quantum 
mechanics [20-27] to solve analytically the Hulthen superpotential partner (/ + l)th member 
for non-zero angular momentum states given by Ref.[18]: 



1-1(1 + 1) 



h z 5 

2mZe 2 



e~ Sr h 2 5 2 e~ 2Sr 

+ ^r*(* + *) „ i = 0,1,2, 



l-e~ Sr 2m v ' (l -e- Sr f 

(2) 



which is for s— states is not shape invariant [28] and leads to the usual Hulthen potential 
(1). On the other hand, Eq.(2) can be rearranged as 

t r(H) / \ „ „ 2r e- Sr l(l + l)h 2 5 2 e~ 5r 

VifUr) = V (l+1) (r) = -Ze 2 5 j- + - — y, 3 

eff\ ) ) l- e ~ 5r 2m (i_ e -^) 2 v; 

which is well-known as the approximate Hulthen effective potential introduced by Greene 
and Aldrich [29]. 1 It is of much concern to see that for small values of 5, Eq.(3) becomes 
the effective Coulomb potential given by 

Or,*»0)^, = -^ + «£i£. (4) 

In the case of Coulomb potential, the Hamiltonian corresponds to the addition of an appro- 
priate barrier potential and the so-called degenarcy is recovered as a natural consequence 
[18]. At small values of r, the Hulthen potential, behaves like a Coulomb potential whereas 
for large values of r it decreases exponentially so its capacity for bound state is smaller than 
the Coulomb potential. The coulomb problem is analytically solvable for all energies and 
angular momenta. Because of the similarity and points of contrast mentioned above between 
Coulomb and Hulthen potentials, it may be of much interest to use the Hulthen superpo- 
tential partner, Eq.(2), to generate their eigenvalues and eigenfunctions in the framework 
of the PT— symmetric quantum mechanics by employing Nikiforov-Uvarov (NU) method. 
The outline of the paper is as follows: In section II, we solve the Schrodinger equation SE 
with Hulthen superpotential for its energy eigenvalues and eigenfunctions. In section III, we 
consider the / = and I ^ cases and compare with the other works and methods. Finally, 
section IV is for our conclusions. 



: In Ref.[30] the Hulthen effective potential is taken as V { e f f (r) = Sj^w + ^^^f in 
atomic units. 
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II. POLYNOMIAL SOLUTION OF THE HULTHEN POTENTIAL 



The NU method provides us an exact solution of non-relativistic SE for certain kind of 
potentials [19]. The method is based upon the solutions of general second order linear 
differential equation with special orthogonal functions [31]. For a given real or complex 
potentials, the SE in one dimension, which is a PT— symmetric equation, is reduced to 
a generalized equation of hypergeo metric type with an appropriate s = s(x) coordinate 
transformation. Thus, it takes the form: 

r(s) + ^'(s) + 4r-Ms) = o, (5) 

cys) a \ s ) 

where a(s) and a(s) are polynomials, at most of second- degree, and f(s) is of a first-degree 
polynomial. To find a particular solution for SE by separation of variables, we use the 
transformation given by 

= <Ks)y(8). (6) 

This reduces SE into an equation of hypergeometric type: 

a(s)y"(s) + r( S )y'(s) + \y(s)=0, (7) 

where is found to isatisfy the condition (j)' (s) / 4>(s) = n(s)/a(s). Further, y(s) is the 
hypergeometric type function whose polynomial solutions are given by Rodrigues relation 

where B n is a normalizing constant and the weight function p(s) must satisfy the condition 
[19] 

Hs)p(s))' = r(s)p(s). (9) 
The function ir(s) and the parameter A required for this method are defined by 



1 



vw-f M y_ g(8)+fa7(s)j 



and 



A =k + TT , (s). 



(11) 



Here, ir(s) is a polynomial with the parameter s and the determination of k is the essential 
point in the calculation of 7r(s). Thus, for the determination of k, the discriminant under 
the square root is being set equal to zero and the resulting second-order polynomial has to 
be solved for its roots k± 2 . Hence, a new eigenvalue equation for the SE becomes 



\ n + nr'(s) + n(n o % "(s) =0, (n = 0, 1, 2, ...) 



(12) 



where 



t(s) = t(s) + 2tv(s), 

and it will have a negative derivative. 

Now, we follow Ref.[32], by rewritting Eq.(2) in a quite simple form as 



= — Sr 



—8r 



V(r) = -V\ 



1-e 



-dr 



1-e 



-5r 



with 



V l = Ze 2 S 



1-1(1 + 1): 



h 2 5 



2mZe 2 

Therefore, using the separation of variables 

1 



and V 2 = —1(1 + 1). 
2m 



1>(r) = -R(r)Y(e,<i>), 



(13) 



(14) 



(15) 



(16) 



we may write the radial part of SE for all angular momentum states as 



2m dr 2 \ 2mr 2 , 

On the other hand, the one-dimensional counterpart of Eq.(17) can be written, in a PT 
symmetric form, as 



(17) 



R"(x) + 



E + 



-8x 



V 2 e 



-2&x 



l- e -sx (i_ e -s x y_ 
5 



R(x) = 0, 



(18) 



and with the assi§gnment s = e Sx , then it becomes 



d 2 R{s) 1 dR{s) 2m 



ds 2 s ds Ti 2 5 2 s 2 



E + 



V lS V 2 s 2 



i-s a-s) 



R(s) = 0, 



(19) 



and also introducing the given dimensionless parameters 



2mE . . n 2mV 1 ,„ . 2mV 2 , 

^->0 (E<0), (3=-^ ((3>0), l = ~ri (7>0), 



n 2 5 2 v n " n 2 P ^ n 1 h 2 5 2 

finally it leads into the following simple hypergeometric form given by 



(20) 



d 2 R{s) (1 - s) dRjs) ^ 1 



x 



- (e + (3 + 7) s 2 + (2e + (3) s - e R(s) = 0. (21) 



ds 2 ' s(l-s) ds ' [s(l-s)] 2 
Hence, comparing the last equation with the generalized hypergeometric type, Eq.(5), we 
obtain the associated polynomials as 

?(s) = l-s, a(s) = s(l-s), a(s) = -(e + (3 + 1 )s 2 + (2e + (3)s-e. (22) 

When these polynomials are substituted into Eq.(10), with cr'(s) = 1 — 2s, we obtain 

tt(s) = ~ ± -y^(l + 4e + 4/3 + 4 7 + 4fc) s 2 - 4 (/? + 2e + k) s + 4e. (23) 



Further, the discriminant of the upper expression under the square root has to be set equal 
to zero. Therefore, it becomes 



A = [4 {(3 + 2e + A;)] 2 - 4 x 4e (1 + 4e + 4/3 + 4 7 + Ik) = 0. 



(24) 



Solving Eq.(24) for the constant k, we get the double roots as A; +i _ = —(3 ± \Je (1 + 47), and 
substituting these values for each into Eq.(23), we obtain 



7T S = ± - < 

y ' 2 2^ 



[(2^ - v/TT47) s - 2^] ; for k + = -(3 + ^/e(l + 4 7 ), 

[(2^+ V ^4^)s-2 V ^] ; for A;_ = -(3 - (1 + 4 7 ). 

Making the following choice for the polynomial ir(s) as 

s 1 



(25) 



7T(s) = 



2 2 



(2^ + ^/l + 4 7 ) s - 2^ 



(26) 



gives the function: 
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t(s) = 1 -2s- ^(Ve+ VI + 47J «-2v^J , (27) 

which has a negative derivative of the form r(s) = — (2 + 2^/e + ^/l + 4^yj . Thus, from 
Eq.(ll) and Eq.(12), we find 



and 



A = -f3 - l - (l + 2y/l) (l + ^1 + 47) , 
A n = - ^2 + 2^1 + ^1 + 47) n - n(n - 1). 



After setting A„ = A and solving for e, we find: 



l + 2n 



(n(n + 1) + (3) 



(28) 



(29) 



(30) 



2 1 + 2n + VI + 47. 
which is exactly Eq.(26) in Ref.[32] for the deformed Woods-Saxon potential if one lets 

q — — 1 and a = 1/5. Therefore, substituting the values of e, (3 and 7 into Eq.(30), one can 

immediately determine the Hulthen's exact energy eigenvalues E n i + i as 



E. 



(H) 
n,l+l 



2x2 



' 2m 



l + 2n (n(n + 1) + ) 



1 + 2n + y^TT 
Therefore, substituting, Eq.(15) into Eq.(31), one gets 

ie 2 Z/ft 2 ) ( n + / + i) 



8my 2 



< n < 00. 



(31) 



^71,1+1 



2m 



< n < 00. 



n + 1 + l 2 - • (32) 

for / + 1 Hulthen superpotential. Following Ref.[32], in atomic units (h — m — c — e — 1) 
and for Z — 1, Eq.(32) turns out to be 



(*0 

n,Z+l 



n + l + 1 



or 



(*0 



(ra + Z + 1) 



n + l 



1 2 



< n < 00, Z = 0,1,2, ... 



n = n+l, / = 0,1,2,.. 



(33) 



(34) 



n + l 2 

which is exactly the same result obtained by other works (cf. e.g., Ref.[33], Eq.(78)) if / 
is set equal to zero. The above equation indicates that we deal with a family of Hulthen 
potentials. 2 Equation (32) agrees with Eq.(5) in Ref.[18] for / 7^ case, and Eq.(ll) in 



2 The critical screening S c , at which E n = 0, is defined, in atomic units, by <5 C = 2/(n + 1 + l) 2 . 
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Ref.[30] for I — case. Of course, it is clear that by imposing appropriate changes in the 
parameters 5 and V±, the index n describes the quantization for the bound energy states. 
In addition, if the parameter V2 in Eq.(31) is adjusted to zero, solution reduces to the form 
obtained for the standard Hulthen potential without a barrier term (cf. e.g., Eqs.(14) and 
(15) with / = case). 

Let us now find the corresponding wavef unctions. Applying the NU method, the polyno- 
mial solutions of the hypergeometric function y(s) depends on the determination of weight 
function p(s) which is found to be 

p(s) = (1 - sy-h 2 ^; 77 = 1 + y/l + 4 7 . (35) 

Substituting into the Rodrigues relation given in Eq.(8), the eigenfunctions are obtained in 
the following form 

Vn, q {s) = C n (l - 8 )-<*-%-*T<£. [(1 - sT^ 1 s n+ ^] , (36) 

where C n stands for the normalization constant and its value is l/n\. The polynomial so- 
lutions of y n (s) are expressed in terms of Jacobi Polynomials, which is one of the classical 
orthogonal polynomials, with weight function (1 — sY'^s 2 ^ in the closed interval [0, 1] , 
yielding A n P^ 2v/ ^'' ? ~ 1 )(i — 2s) [31]. Finally, the other part of the wave function in Eq.(6) is 
found to be 

0( s ) = (1 - s f s v^ ^ = v / 2 . (37) 

Combining the Jacobi polynomials and 0(s) in Eq.(36), the s— wave functions (I = 0) could 
be determined as 

R n (s) = D n s^(l - S yP^-V(i _ 2s ), (38) 
with s = e~ Sx and D n is a new normalization constant. 
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III. CONCLUSIONS 



The exact solutions of the radial SE for the Hulthen potential with the angular mo- 
mentum I = and I 7^ are found by using NU method. Eigenvalues and eigenfunctions 
obtained from the real form of the potential are computed. Therefore, the wave functions 
are physical and energy eigenvalues are in good agreement with the results obtained by the 
other methods. In this regard, Figure 1 shows the variation of the Hulthen potential with 
r for S-, P— , and D— states. Figure 2 plots the variation of the Hulthen potential with 
r for the S— state with various screening parameters 5 = 0.002,0.01 and 0.1. Figures 3 
and 4 show the variation of the energy eigenvalues with respect to the quantum number 
n for S— and P— states with the a chosen values of screening parameter 5 = 0.002, 0.05 
and 0.2, respectively. On the other hand, Table 1 shows the bound energy eigenvalues of 
the Hulthen potential as a function of 5 for various quantum numbers of S'-state. These 
results are compared with other works [11,33]. Further, the bound energy eigenvalues as a 
function of 5 for the states 2p (n — 0, 1 — 1) and 3d (n — 1,1 — 1) [30] are given in Table 2. 
Comparison of our results with numerical data of Refs.[3,30] is also given. However, since 
the form of the potential used in our work has a different form than the potential form used 
by Ref.[30] in Eq. (14), we have to add a perturbation term, AE, to our calculations in 
order to substitute the small differences in [18]. Better results have been obtained for the 
2p state for small values of 5 since the effective potential (3) becomes closer to the original 
Hulthen potential (1) and for small / the contribution of this angular momentum term in 
potential is also small. Therefore, if all the parameters of potential remain purely real, it 
is clear that all bound energies E n with n > represent a negative energy spectrum [32]. 
We also point out that the exact results obtained for the Hulthen potential may have some 
interesting applications in the study of different quantum mechanical systems and atomic 
physics. 
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FIGURES 

FIG. 1. Variation of the Hulthen potential as a function r. The curves are shown for screening 
parameter 5 = 0.2 for the S—,P—, and D— states. 

FIG. 2. Variation of the Hulthen potential as a function r. The curves are shown for S— state 
with various values of the screening parameter S = 0.002, 0.01 and 0.1. 

FIG. 3. The variation of the energy eigenvalues with respect to the quantum number n. The 
curves shown are for S— state with various values of the screening parameters 5 = 0.002, 0.05 and 
0.2. 

FIG. 4. The variation of the energy eigenvalues with respect to the quantum number n. The 
curves shown are for P— state with the same screening parameters as in Figure 3. 
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TABLES 

TABLE I. The S-states energy eigenvalues of the Hulthen potential for several values of screen- 
ing parameter 5. 
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TABLE II. Energy eigenvalues as a function of the screening parameter 5 for the states 2p and 

3d. 



State 


5 


-E n ,i [3] a 


-E n ,i [30] b 


Our work c 


2p 


0.025 


0.112760 


0.1127605 


0.1128125 




0.050 


0.101042 


0.1010425 


0.1012500 




0.075 


0.089845 


0.0898478 


0.0903125 




0.100 
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0.0450000 


3d 
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0.043601 


0.0437069 


0.0437590 




0.050 


0.032748 


0.331645 


0.0333681 




0.075 


0.023010 


0.0239397 


0.0243837 




0.100 


0.014433 


0.0160537 


0.0168056 



Variational method. 
b Numerical integration. 

c The small difference in results is because the potential form used by Ref.[30] for I / part is 
different than our form. 
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